Ferromagnetism in Co 7 (Te03)4Br 6 : A byproduct of complex antiferromagnetic order 

and single-ion anisotropy 
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Pronounced anisotropy of magnetic properties and complex magnetic order of a new oxi-halide 
compound Co7(Te03)4Br6 has been investigated by powder and single crystal neutron diffraction, 
magnetization and ac susceptibility techniques. Anisotropy of susceptibility extends far into the 
paramagnetic temperature range. A principal source of anisotropy are anisotropic properties of 
the involved octahedrally coordinated single Co 2+ ions, as confirmed by angular-overlap-model cal- 
culations presented in this work. Incommensurate antiferromagnetic order sets in at Tjv=34 K. 
Propagation vector is strongly temperature dependent reaching ki = (0. 9458(6), 0, 0.6026(5)) at 30 
K. A transition to a ferrimagnetic structure with k2=0 takes place at Tc=27 K. Magnetically or- 
dered phase is characterized by very unusual anisotropy as well: while M — H scans along b-axis 
reveals spectacularly rectangular but otherwise standard ferromagnetic hysteresis loops, M — H 
studies along other two principal axes are perfectly reversible, revealing very sharp spin flop (or spin 
flip) transitions, like in a standard antiferromagnet (or metamagnet). 

Altogether, the observed magnetic phenomenology is interpreted as an evidence of competing mag- 
netic interactions permeating the system, first of all of the single ion anisotropy energy and the 
exchange interactions. Different coordinations of the Co 2+ -ions involved in the low-symmetry C2/c 
structure of Co7(Te03)4Br6 render the exchange-interaction network very complex by itself. Tem- 
perature dependent changes in the magnetic structure, together with an abrupt emergence of a ferro- 
magnetic component, are ascribed to continual spin reorientations described by a multi-component, 
but yet unknown, spin Hamiltonian. 

PACS numbers: 12.34, 56.78 



I. INTRODUCTION 



Enormous diversity of magnetic phenomena relies on 
equal diversity of interactions permeating real magnetic 
materials. In cases of one interaction dominating by 
far over the others (typically valid for exchange inter- 
action) an elementary insight into magnetism of such a 
system can indeed be acquired on basis of a simple, one- 
interaction Hamiltonian. However, any profound knowl- 
edge of magnetism of even such a simple system implies 
taking into account other interactions present in the sys- 
tem, whatever weak they could be. For example, basic 
static properties (like magnetization, susceptibility, mag- 
netic structure...) of long-range ordered ferromagnets 
and antifcrromagnets, especially if they rely on strong 
exchange interaction, can be well-interpreted within one- 
interaction models. Understanding dynamic features of 
the same systems, like their spin excitation spectra for 
example, requires however at least magnetocrystalline 
anisotropy (thus, in turn, several inevitably involved spe- 
cific interactions) to be included as well. The most in- 
teresting situation, for fundamental research as well for 
applications, arises in cases involving presence of mag- 
netic interactions competing each other in size and/or 
in sign while residing on energy scale of thermal excita- 
tions within the usual experimental window 2 — 300 K. 
Antiferromagnetic exchange competing with single ion 



anisotropy can, for example, render the certainty of long- 
range ordering questionable, as shown decades ago for 
archetypal antiferromagnet^ N1F2. In a more recently 
developed general framework of Quantum Magnetism 2 
the subject of competing interactions is recognized as 
a key ingredient introducing the quantum phase tran- 
sition point into the respective phase diagrams. In an- 
other aspect, competing interactions are responsible for 
complex/incommensurate magnetic structures and emer- 
gence of ferromagnetism in the ground states of many 
nominally antiferromagnetic systems^. In this category 
there is a particularly interesting group of systems reveal- 
ing spin reorientations taking place within the system- 
specific temperature intervals. The prominent exam- 
ples are the rare-earth elements Dysprosium (Dy) and 
Terbium (Tb)^>, transition-metals sesquioxides^ (notably, 
hematite- 6 , cv-Fe203) and orthoferritesi From applica- 
tive side there is renewed interest for these systems, par- 
ticularly for orthoferritesi^: spins in these systems are 
subject to combined effect of antiferromagnetic exchange 
and magnetocrystalline anisotropy which enables ultra- 
fast spin manipulation 12 . Ultrafast dynamics is a key 
issue for exchange-bias devices. 

In this article we present magnetic structure and 
properties of a recently discovered^ magnetic system, 
Coy(Te03)4Br6, which, as shown herewith, reveals re- 
markable manifestations of competing magnetic interac- 
tions. In its ground state there is a complex noncolin- 
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ear long-range magnetic order while pronounced mag- 
netic anisotropy characterizes both paramagnetic and 
magnetically ordered phase. In the temperature range 
27K — the incipient incommensurate magnetic or- 
der continuously changes by cooling, ending up at 27K 
by an abruptly intruding ferromagnetic component along 
one crystallographic axis. In the spin reorientation, un- 
derlying emergence of ferromagnetism, antiferromagnetic 
backbone remains conserved: we show that the key hall- 
mark of antiferromagnetic order, spin flop/flip transi- 
tion, keeps characterizing magnetism of Co7(Te03)4Br 6 
in measurements along other two principal axes. The re- 
sults are interpreted by invoking competition of exchange 
interactions primarily with single ion anisotropy energy 
but also with several other possible sources of magne- 
tocrystalline anisotropy. 

The article is organized as follows. The involved exper- 
imental techniques are introduced in Section [TTJ In Sec- 
tion lllll results of neutron diffraction, dc/ac susceptibility 
as well as of magnetic hysteresis/M-H studies are pre- 
sented. Results of theoretical modeling of paramagnetic 
susceptibility and of the ground state properties, based 
on angular-overlap-model calculations, are also presented 
in Section IIIII Both experimental and theoretical results 
provide firm evidence of exchange interactions and sin- 
gle ion anisotropy energy ruling magnetic properties of 
Co7(Te03)4Br6. In Scction HVl these results are discussed 
in a more general framework of competing interactions of 
various sorts, knowledge of which was accumulating dur- 
ing the decades. In Section [V] appropriate conclusions 
have been presented. 



II. EXPERIMENTAL DETAILS 

The single crystals of Co7(Te03)4Br6 were synthesized 
via chemical vapor transport reactions. The details of the 
synthesis can be found elsewhere^. The single crystals 
grow in platelet geometry and the typical samples used 
in the magnetization/susceptibility studies had approxi- 
mate dimensions 4.0 x 2.0 x 0.1 mm 3 . The plane of the 
platelet-like samples corresponds to the crystallographic 
&c-plane. 

Neutron powder diffraction data have been collected 
in the temperature range 3 K - 45 K on the DMC instru- 
ment at SINQ, Paul Scherrer Institute, Villigen, Switzer- 
land, with neutron wavelength of 2.453 A and 4.2 A. In- 
commensurate (ICM) wave vector has been determined 
on the single crystal diffractometer TriCS at SINQ with 
neutron wavelength of 2.32 A using the area detector. 

Magnetization/susceptibility measurements were per- 
formed on oriented single crystals in applied magnetic 
fields directed parallel to the crystallographic a*-, b-, 
and c-axis. Samples were oriented by the use of X-ray 
diffractometer. The choice of a*-axis sample orientation, 
instead of the preferable a-axis one, is imposed by sample 
morphology. For DC magnetization studies up to 5.5 T a 
Quantum Design superconducting quantum interference 



device (SQUID) magnetometer was used, covering the 
temperature range 2 - 300 K. AC susceptibility studies 
were performed using a CryoBIND ac susceptibility sys- 
tem employing measuring ac field of 3 Oe and frequency 
of 430 Hz. 



III. RESULTS 
A. Neutron diffraction 

Refinement of the neutron powder pattern collected at 
45 K confirmed that the sample is a single phase with 
a monoclinic C2/c space group and unit cell parameters 
a = 20.590(2) A, b = 8.4998(7) A, c = 14.631(1) A, (3 = 
125.202(6)°, in close agreement with the original struc- 
tural workS. As pointed out therein the crystal structure 
can be described as layered in the frc-plane. The layers 
are built of networks comprising three types of Co 2+ - 
ion distorted octahedra ([Co(l)04Br 2 ], [Co(2)04Br 2 ], 
[Co(3)04Br2]) and [TeOsE] tetrahedra, while the layers 
are interconnected along the a-axis by the fourth Co 2+ - 
ion octahedron type, [Co(4)02Br 4 ]. 

The collected neutron powder diffraction patterns con- 
firmed the existence of two sequential magnetic orderings, 
as claimed by the original work on basis of magnetic sus- 
ceptibility studies^. Below Tn — 34 K weak magnetic 
peaks occur. They cannot be indexed as simple multi- 
ples of the crystal unit cell revealing an incommensurate 
wave vector. Moreover, the incommensurability appears 
in two directions, a* and c*, which has been clarified by 
a single crystal study using the 2D detector of TriCS. 
The position of the ICM peaks is strongly temperature 
dependent and in Fig. [T] we show the difference neutron 
powder diffraction patterns T-45 K in the range T= 25 - 
32 K. Each pattern presented in Fig.[T]has been analyzed 
using the profile match option of the Fullprof program 9 
and the two ICM components (k x , 0, k z ) of the wave vec- 
tor has been refined. The temperature dependence of k x 
and k z is shown in Fig. [21 

At 27.5 K a second set of additional strong magnetic 
peaks appears, coexisting in a short temperature interval 
with the ICM peaks. Further temperature lowering leads 
to weakening of the ICM peaks and their transformation 
into diffuse scattering, as shown in the bottom pattern 
of Fig. [TJ Below 26 K only the second set of magnetic 
peaks remains. 

This set can be indexed with the wave vector k 2 = 0. 
Some of the new magnetic peaks (i.e. (200), (110)) over- 
lap with the nuclear reflections implying ferromagnetic 
contribution, others (i.e. (001), (-201)) appear at the po- 
sitions extinct in the paramagnetic pattern, as expected 
from an antiferromagnetic component. The systematic 
extinctions observed in the 5 K - 45 K magnetic differ- 
ence pattern reveal that the C-translation and the glide 
plane c not combined with the time reversal are retained 
in magnetic symmetry {hkl : h+k = 2n and hOl : I ^ 2n). 

Representation analysis implemented in the Fullprof 
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FIG. 1: (Color online) Difference neutron powder diffraction 
T - 45 K of Co 7 (Te0 3 )4Br6 with T in the range 25 - 32 K 
(DMC instrument, A = 4.2 A). 




T[K] 



FIG. 2: (Color online) Temperature dependence of the ICM 
magnetic vector components refined from the DMC data, A = 
2.453 A. 



program^ has been used to determine the k2 = mag- 
netic structure. The Fourier coefficients describing pos- 
sible spin configurations can be written as linear com- 
binations of irreducible representations (IR) of the wave 
vector group (little group). The magnetic representations 
for the 4a and 8f sites, occupied by Co(4) and Co(l-3), 
can be decomposed in IR's: 

r(4a)= 3 ri + 3 r 3 

r(8f)= 3 ri + 3 r 2 + 3 r 3 + 3 r 4 

Note that there are four ions of the Co(4) set and eight 
ions in each of the Co(l), Co(2) and Co(3) sets. In the T\ 
and T 3 representations magnetic moments of the Co(l)- 
Co(4) sets may attain independent values and directions, 
while the moments of the ions within the same site are 
constrained by the symmetry relations presented in Ta- 
ble |TJ So all together there are twelve independent pa- 
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FIG. 3: (Color online) Observed 5 K - 45 K magnetic 
difference pattern, calculated and difference patterns of 
Co7(Te03)4Brs denoted by crosses, red solid, and green dot- 
ted lines, respectively. Inset: The choice of the XY Z — a* be 
orthogonal system and morphology of the single crystal used. 



TABLE I: Irreducible representations of the wave vector 
group for \i2~ in the space group C2/c. The notation Co(ij) 
is used with the index i=l,...,4 labeling the set, and the index 
j=l,...,4 labeling the ions within the site. The ions of the 
same site have coordinates: j=l (x y z), j=2 (1-x y -z+1/2), 
j=3 (-x+1/2 -y+1/2 1-z), j=4 (x 1-y z+1/2). The ions gen- 
erated by the C-lattice translation have the same magnetic 
moment as the generating atom and are therefore omitted. 
The u,v,w coefficients are fixed by the symmetry for one set, 
but are independent for different sets. 



ion r a r 2 r 3 r 4 

4 a 

Co(41) uvw - uvw 

Co(43) -u v -w - u -v w 

i=l,3 8 f 

Co(il) uvw uvw uvw uvw 

Co(i2) -u v -w -u v -w u -v w u -v w 

Co(i3) uvw -u -v -w uvw -u -v -w 

Co(i4) -u v -w u -v w u -v w -u v -w 



rameters (m x , rn y and m z of four sets of Co 2+ ions) in 
Ti and r 3 . In T2, T4 IR's the magnetic moments of the 
Co(4) set must be zero. 

Refinement of the models has been performed with 
Fullprof^. The best agreement with experimental data 
(Rm — 6.6%, Fig. [3]) is obtained for a 3-dimensional 
canted Ti model presented in Table [Til an d Fig. [4j The 
a* be orthogonal coordinate system defined in Inset to 
Fig. [3] has been used. The Co moments reach the values 
of 4.21(7) /VCo(l), 4.4(1) Ms/Go(2), 3.8(1) aVCo(3) 
and 4.5(1) ^s/Co(4) at 5 K. These values are larger than 
the spin only component (3 /is) of the Co 2+ ion confirm- 
ing incomplete quenching of the orbital moment for all 
cobalt ions in this compound 8 . 
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TABLE II: Refined parameters for the \i.2~ magnetic struc- 
ture. M[/Xb] is the ordered magnetic moment, with the M. xyz 
components defined in the a* be orthogonal coordinate system. 
The Co(l-4) sets are represented by the ions j—1 (x y z) and 
j=2 (1-x y -z+1/2). a [deg] is the canting angle between the 
magnetic moments of the 1st and the following ions. 



Ions 


M. 


M, 


M z 


Ms /Co 


aii 


Co(ll) 


2.64(9) 


1.86(6) 


2.77( 9) 


4.25(7) 





Co(12) 


-2.64(9) 


1.86(6) 


-2.77 { 9) 


4.25(7) 


125(2) 


Co(21) 


-3.08(7) 


-1.1( 1) 


3.08(7) 


4.4(1) 


95(2) 


Co(22) 


3.08(7) 


-1.1( 1) 


-3.08(7) 


4.4(1) 


97(2) 


Co(31) 


-0.43(8) 


-3.1( 1) 


2.22(8) 


3.8(1) 


93(2) 


Co(32) 


0.43(8) 


-3.1( 1) 


-2.22(8) 


3.8(1) 


131(2) 


Co(41) 


-3.51(8) 


1.33( 9) 


2.48(8) 


4.5(1) 


90(2) 


Co(42) 


3.51(8) 


1.33( 9) 


-2.48(8) 


4.5(1) 


75(2) 




FIG. 4: (Color online) The ac-projection of the low- 
temperature (5 K) magnetic structure of Co7(TeOa)4Br6. 
The four crystallographic sets of Co 2+ ions are shown by vi- 
olet - Co(l), red - Co(2), green - Co(3) and orange - Co(4). 
The Co(21) and Co(32) ions, as well as Co(22) and Co(32), 
superimpose on this projection. The pairs Co(ij)/Co(ij+l) 
have opposite and M z magnetic components, but the same 
My component. The sign of M y is shown near the symbol of 
each ion. Spatial orientations of local magnetic moments are 
illustrated by three-dimensional vectors. 



The Co7(Te03)4Br 6 system could be identified as a 
canted antiferromagnet (weak ferromagnet) - the angle 
between moments of the similar magnitude is smaller 
than 180°, giving rise to a ferromagnetic component. 
Each Co set (i — 1,...,4) has a different magnitude of 
the My component. Ferromagnetic components associ- 
ated with the Co(l) and Co(4) sets point in one direc- 
tion, while those associated with Co(2) and Co(3) point 
in the opposite direction. The net moment along the 
b axis given by the sum of the M y values in Table HI1 
amounts to « 0.5/is/Co, which is in close agreement 
with the remanent moment obtained from the magneti- 
zation measurements (see, section D). As elaborated in 
the discussion section, the complex non-collinear canted 
magnetic structure of Co7(Te03)4Br6 (Fig.QJ stems from 
different involved interactions, first of all from the com- 
peting single ion anisotropy and exchange interactions. 
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FIG. 5: (Color online) Main Panel: Temperature dependence 
of real component of ac susceptibility (3 Oe, 430 Hz) for mea- 
suring field oriented along three principal axes. DC suscepti- 
bility in small applied field (100 Oe) provides almost identical 
result. Black curve represents the Curie plot for free S = 3/2, 
g=2.5, ions (Curie constant C=20.5 emuK/mol) Inset Top: 
Product xdc ■ T vs. temperature. Inset Bottom: Susceptibil- 
ity inverse x~dc vs - temperature. 



B. Magnetic susceptibility 

1. High-temperature dependence 

Results of ac and dc susceptibility measurements in 
three crystallographic directions a* , b and c are shown in 
Fig-O Extending the original susceptibility report^ these 
results document pronounced susceptibility anisotropy 
characterizing Co7(Te03)4Brg in a broad temperature 
range. For all three sample orientations the results reveal 
the presence of magnetic transition at Tn = 34 K, while 
the measurement along one axis only (6 axis) documents 
the presence of an additional, very pronounced transition 
at Tq — 27 K, in agreement with the original study. The 
present study shows that in the paramagnetic range the 
particular directional susceptibilities are remarkably dif- 
ferent, the difference extending far above the transition 
temperature range. Unlike very anisotropic susceptibil- 
ity the effective <?-factor value, as determined from the 
high-temperature limit of the Curie- Weiss (CW) plot, 
was found to be pretty isotropic in the room temperature 
range. The values for g were found to be very close to 
g = 2.5 (with S=3/2) for all three sample-to-field orienta- 
tions, in full accordance with the original powder-sample 
data-. The value of the Weiss parameter 9 of the CW plot 
varies depending on the chosen axis. We elaborate be- 
low that the behavior of the susceptibility is dominantly 
influenced by the single-ion anisotropy so the values of 9 
cannot be used to estimate the type and the strength of 
the involved exchange interactions. From the experimen- 
tal side we note that there were only marginal sample- 
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FIG. 6: (Color online) Temperature dependence of out-of- 
phase (imaginary) component of ac susceptibility for three 
principal axes. Above the transition range there is over- 
lap of the data points for all three measurement directions. 
The data for x (Fig- © aQ d x" were taken within the same 
experimental runs. Inset: Zoomed transition temperature 
range. Residual signal in measurement along c-axis (green 
symbols) is ascribed to imperfect sample alignment and/or 
non-vanishing cross-talk between the in- and out-of-phase 
components in the phase sensitive detection. 



to-sample and batch-to-batch variations: The reported 
results thus rely only on the intrinsic crystal structure of 
the compound. 

Temperature dependence of the imaginary component 
of susceptibility is shown in Fig. [5] These data are strik- 
ing in two aspects. Firstly, imaginary susceptibility sig- 
nal is present in the direction of the 6-axis only. Secondly, 
the size and the sharpness of the imaginary susceptibil- 
ity peak is, to the best of our knowledge, not observed 
in other magnetically ordered systems: under particular 
conditions (see Fig. [6]) imaginary peak is more than two 
times bigger than the peak in real susceptibility. Tra- 
ditionally, peak in imaginary susceptibility is ascribed 
to dissipative ferromagnetic-domain dynamics, setting in 
simultaneously with the formation of the domains im- 
mediately below T c . Due to its extraordinary properties 
imaginary susceptibility of Co7(Te03)4Brg will be sub- 
ject of a separate publication. 



2. Susceptibility in magnetically ordered phase 

At Tjv = 34 K a long-range magnetic order sets in. 
Both magnetic susceptibility and neutron diffraction data 
are consistent with antiferromagnetic ordering. Rela- 
tive to the susceptibility maxima at Tjv there is a rapid 
drop of both c— and a*— axis susceptibility by cooling 
below Tjv (Fig. [5]), while 6— axis susceptibility simul- 
taneously builds up. These findings are incompatible 
with simple uniaxial magnetic AF, which obviously does 
not take place in Co7(Te03)4Bi'6. Three-dimensionally- 
canted and incommensurate order, documented by neu- 



tron diffraction results presented above, is fully compat- 
ible with susceptibility data. One notes that at Tc, 
marked by a very sharp peak of b— axis susceptibility, 
c— and a*— axis susceptibility exhibit a sudden drop 
(Fig. [5]) . The most natural explanation is that at Tc an 
abrupt magnetic moment reorientation takes place such 
that a growing ferromagnetic component builds up at the 
expense of magnetic moments participating in the ICM 
ordering at higher temperatures. In accordance with the 
low-temperature magnetic structure (Fig. 2]) the latter 
observation shows that in magnetic ordering all magnetic 
degrees of freedom participate cooperatively, i.e., even if 
there exist separate magnetic fractions they contribute 
in mutual accord into the total susceptibility. One has 
to point out that low-field susceptibility studies reveal 
a remarkable feature of the system that one could not 
figure out from the crystallographic or low-temperature 
magnetic structures alone: in spite of complexity of mag- 
netic structure susceptibility data shows that there is a 
distinct axis, aligned approximately along c— axis, which 
defines the direction of preferable spin orientation^. (A 
more precise determination of the preferred axis orienta- 
tion is presented in Subsection IIII Dl ) A fact that in the 
low-temperature magnetic structure the spins are aligned 
along very different directions, thus not along c- nor any 
other axis, is a consequence of different competing inter- 
actions ruling the spin geometry in the ground state. 



C. Anisotropic magnetic properties in Angular 
Overlap Model calculations 

Angular Overlap Model (AOM) calculations were per- 
formed to estimate the single ion magnetic anisotropy 
arising from the crystallographically inequivalcnt 
cobalt (II) centres. AOM parameters for the Co-0 
and Co-Br bonding interactions were estimated from 
values of the ligand field splitting parameters tabulated 
for octahedral Co(H20)g + and tetrahedral CoBr^ - 
complexe o 30 i 31 assuming e w = 0.2e CT . The parameter 
e CT was assumed to vary with distance 3 ^ as a function 
of 1/r 5 and e n as a function of 1/r 6 . The Racah 
and spin-orbit coupling parameters were fixed at 80 
% of their free-ion values, and the orbital Zeeman 
interaction reduced accordingly. The AOM matrices 
were constructed using LIGFIELD 29 with the AOM 
parameters and angular coordinates calculated from 
structural data as input, employing all 120 functions of 
the 3d 7 electronic configuration. For the calculation of 
magnetic susceptibility curves, the AOM matrices were 
imported into the program MagProp 33 . 

The magnetic moment per ion was calculated from the 
expression 

M E.(-frM-fr) 
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and the paramagnetic molar susceptibility from 



X P = Na : 



A4_ 
B 



(2) 



In these equations B designates the external magnetic 
field, ks the Bolzmann constant and Na Avagadros num- 
ber. The sum is over the n eigenstates of the Hamilto- 
nian whose energies are designated by the symbol E n . 
The derivative in equation (1) was found according to 
the Hellman-Feynman theorem, 



dEn 

dB 



dB 1 Vnl 



(3) 



The onset of short-range ferromagnetic order was in- 
corporated by expressing the susceptibility as 



XM X P 



(4) 



where A is the molecular field parameter. 

It is seen from FigJ7]that both experiment and AOM- 
based susceptibility calculations identify the component 
of the susceptibility tensor along the c-axis to be larger 
than the components along the a*- and b-axes. This 
would suggest that single-ion anisotropy indeed plays 
a major role in determining the observed magnetic 
anisotropy. Note, however, that the calculated relative 
magnitudes of the susceptibility tensor along the a*- and 
b-axes are not in accordance with experiment. This could 
arise from our rather crude estimate of the AOM bonding 
parameters or may reflect the fact that no magnetic in- 
teraction paths have been introduced into the model. In 
upgrading the model with appropriate interaction paths 
it would be natural to assume that the intralayer bc-plane 
exchange coupling dominates over the coupling along the 
out-of-plane a*-axis direction. Such corrections would 
certainly make the results of the model calculations in 
closer agreement with experiment. 

As a result of the low-symmetry ligand field and spin- 
orbit coupling the 4 Ti g (0h,) ground term is split into 6 
Kramers doublets that are well-separated in energy. In 
particular the first-excited Kramers doublet lies between 
100 and 300 cm" 1 above the ground state for the four 
crystallographically inequivalent cobalt centres. At suffi- 
ciently low temperatures, therefore, the electronic struc- 
ture may be approximated as a pseudo- S=l/2 system, 
this being an approximation that is commonly employed 
for octahcdrally co-ordinated cobalt (II) complexes^. 

The g 2 tensor was calculated according to a method 
described in detail previously^ in which the energies of 
the states of the lowest lying Kramers doublet are mod- 
eled by the eigenvalues of the S=l/2 spin-Hamiltonian, 



H s = 0B ■ g • S. 



(5) 



The orientation of the g 2 -tensor in the a*bc reference 
frame and the ligand arrangement of four Co 2+ sets are 
shown in Fig [51 
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FIG. 7: Calculated products \ ' T for each of the four Co(i) 
octahedra employing full single-ion Hamiltonian and AOM. 
Calculated results for real single crystal as a whole is shown 
at bottom (plot at left). Results corrected for weak ferromag- 
netic interactions (A = 2, see text) is also shown, (plot at 
bottom right), to be compared with the experimental results 
for X ■ T, Fig[5] 



TABLE III: Calculated components of the g-matrix in the 
a* be reference frame and g' in the eigen coordinate frame. 



9ij 


Co(l) 


3=1, 


,3 


Co(2) 




.,3 


i=l 


5.2251 


1.0467 


2.4724 


2.8165 


0.1922 


-0.1903 


2 


1.0468 


2.1070 


0.5217 


0.1922 


3.4392 


1.8252 


3 


2.4724 


0.5217 


4.1107 


-0.1903 


1.8252 


5.6682 


9jj 


2.2740 


1.7214 


7.4476 


2.9437 


2.2866 


6.6936 


m 


Co(3) 


i=i, 


.,3 


Co(4) 


j=b 


.,3 


i=i 


1.9333 


-1.0734 


0.8493 


3.7126 


0.8268 


-7.2514 


2 


-1.0734 


4.8725 


-2.2980 


0.1097 


0.0012 


-1.2788 


3 


0.8493 


-2.2980 


4.2329 


-1.8364 


-0.9816 


3.9631 


9jj 


2.2383 


1.5742 


7.2262 


0.4946 


-0.4940 


9.3942 



A sizable c-axis susceptibility growth over the Curie 
margin (Fig. [5]) is interpreted therefore as a combined 
effect of single-ion anisotropy and short range ferromag- 
netic correlations. On one side, it would be reasonable 
to assume that these correlations rely on weak ferromag- 
netism that, in principle, could accompany growth of 
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FIG. 8: The alignment of g 2 -tensor (violet ellipsoid) and the 
ligands (oxygen in red, bromine in green) of four different Co 
set in the XYZ—a*bc orthogonal coordinate system. 



the short-range antiferromagnetic order as the temper- 
ature approaches TV from above. This possibility has 
to be abandoned, however, as there is no trace of any 
ferromagnetic order (in the form of magnetic hysteresis, 
peak in imaginary susceptibility, see. Fig. [5]) setting in 
at, or immediately below, TV- Ferromagnetic order, set- 
ting in below Tc, can only partially be related to high- 
temperature ferromagnetic correlations. As shown be- 
low, low-temperature ferromagnetism is manifested along 
b-axis only (thus not along the preferred c-axis) and its 
evolution is primarily correlated with the transforma- 
tion/reorientation of the preformed ICM matrix. 



D. M-H studies 

The most common experimental hallmarks of mag- 
netically ordered ferromagnets and antiferromagnets are 
hysteris and spin-flop transition, respectively, known to 
characterize the respective magnetization vs. field (M-H) 
characteristics. Here we present the results of compre- 
hensive M-H studies on Coy(Te03)4Br6 demonstrating 
presence of both of the mentioned hallmarks: which one 
of the two gets activated depends on the chosen orienta- 
tion of principal axes with respect to the applied dc field. 
In brief, the field component along preferred axis, as im- 
posed by single ion anistropy energy (approximately o 
axis) , activates antiferromagnetic response while the field 
component along 6-axis activates ferromagnetic response. 



1.0 




-1.5 - 



-6 -4 -2 2 4 6 

H„H (T) 

FIG. 9: (Color online) Magnetization (per one Co ion) vs. 
field characteristics for c— axis direction in a broad tempera- 
ture range. M(H) curves reveal almost switching but perfectly 
reversible behavior of magnetization in magnetically ordered 
phase. 



1. Antiferromagnetic response 

In Fig. [9] we first show the M — H scans for the dc field 
applied parallel to the effective preferred axis, c-axis. For 
T > TV the M — H curves are closed (i.e., do not show 
up hysteretic loops) and reveal no saturation for high 
fields indicating a paramagnetic state of the system. For 
Tc < T < TV a sharp magnetization jump shows up 
around 20kOe. Magnitude of the magnetization jumps 
increases as T approaches Tc- The field value at which a 
jump occurs, Hg F , does not change within this temper- 
ature interval. For T < Tc, the magnetization transition 
additionally sharpen and Hg F starts to grow by cooling 
and reaching Hg F —40 kOe at 5 K. Magnetization jumps 
are naturally ascribed to the spin-flop type of spin reori- 
entation. The transition region is itself very narrow - at 
20 K the transition is about 200 Oe wide. 

With the field applied along a*-axis (Fig. [T0|) a very 
similar behavior to the case of H \\ c has been obtained. 
In Fig.[TU]we show just several characteristic M—H scans 
(out of many measured) taken at temperatures below and 
immediately above T c . 

Qualitatively, the M — H characteristics for TT 
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c, a* can be analyzed as a switching/spin reorientation 
(spin-flop) phenomenon superimposed over some field- 
dependent background. The background is attributed to 
coherent rotation of sublatices' magnetization and, be- 
low T c , to contributions from the b- axis ferromagnetism 
(s., next section). The background is more pronounced 
in the case H \\ a* and one notes a drastic change of the 
background slope by cooling below T c , Fig. [TUJ A rather 
pronounced background present in H || a* orientation is 
a consequence of imperfect crystal alignment allowing a 
mixing-in of ferromagnetic H \\ b component. Otherwise, 
the major difference between the M — H characteristics 
measured along c- and a*- axis are about factor of 3-4 
bigger spin- flop fields Hg F , in comparison with H SF at 
the same temperatures. The difference is attributed to 
geometrical reasons: in spin-flop transition the effective 
field component is only the one which is aligned along the 
preferred axisi&. From the spin-flop field values measured 

at 20 K (H§p = 4.7 T and H C SF = 1.3 T) and crystal axes 
geometry (Inset to Fig. [3J one easily determines that the 
axis compatible with the minimum spin-flop field closes 



in the (a*,c)-plane the angle <j> (tan( 



H 



with the 



c-axis. From the latter observation one concludes that 
the effective preffered axis is actually declined from the 
c-axis for the angle <^>(s=s 15°). 

The observed sharpness of the spin-flop transition 
(Fig. [9J represents a direct consequence of a sizable mag- 
netic anisotropy characterizing the system. Generally, 
in a spin flop transition the spins of the sublatices ro- 
tate at Hsf to the direction perpendicular to the field 
(and the preferred axis) direction. By subsequent field 
increase spins are coherently rotated to become aligned 
with the field only at the field value H c . At T = OK, 
the ratio of the two critical fields is expected to obey the 



relatio: 



11 HefO) 



l) 1 / 2 . Here, H E and H A rep- 



Hsf(O) 

resents the mean-field exchange field and the anisotropy 
field, respectively. Hence, in cases with inherently big 
anisotropy, such that Ha approaches He , H c is not much 
bigger than H$f implying almost direct reorientation of 
a sublatice spin into the direction of preferred (i.e., mag- 
netic field) axis. Such a transition is usually referred to 
as spin flip and represents a generic feature of metamag- 
netsii. Co7(Te03)4Br6 may therefore be considered as a 
typical metamagnet system. 

Temperature dependence of the field Hs f is shown in 
Fig. [TTJ Although there is no known generic analytic 
form of Hsf(T), for metamagnets and antiferromagnets 
in general, the observed approximately linear tempera- 
ture dependence is not very commonii. We note how- 
ever that in the case of Coy(Te03)4Br6 spin flip tran- 
sition takes place within the magnetic structure which 
involves a ferromagnetically ordered component (along fa- 
axis, next section) thus deviation from behavior known 
for other antiferromagnets/metamagnets should not be 
any surprising. Also, we point out that, in view of com- 
plex three-dimensional magnetic structure (Fig. 0]) the 
term 'spin flip' cannot be applied literary in its text- 
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FIG. 10: (Color online) Magnetization (per one Co ion) vs. 
field curves for a* direction at three characteristic tempera- 
tures. Spin-flop transitions are marked with arrows. 
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FIG. 11: (Color online) Temperature dependence of spin-flip 
field Hg F , high field saturation magnetization M sat and the 
slope in MH curve at high fields xhf, all for magnetic field 
applied along c-axis. 



book meaning, i.e., to mimic complete reversal of pairs 
originally anti-parallel oriented spins into the direction 
of preferred axis. From the magnitude of magnetization 
jump one cannot associate spin-flipping to any particu- 
lar ion pair belonging to the two sublatices, which could 
perform, hypothetically, spin reversal isolated from the 
rest of structure. Instead, it's more probable that in the 
energy landscape of magnetic structure as a whole there 
are two neighboring minima becoming equal in energy by 
the application of magnetic field Hsf along the effective 
preferred axis. 



2. Ferromagnetic response 

When the magnetic field is applied along the b axis, 
along which the ferromagnetic component exists, remark- 
ably different behavior is observed, Fig. [TJ] Below Tjy a 




-5 -4 -3 -2 -1 1 2 3 4 5 



H H ( T ) 

FIG. 12: (Color online) Main Panel: Magnetic hysteresis 
loops for the fe-axis direction at different temperatures. Mag- 
netization is scaled to one Co ion. Inset: Virgin hysteresis 
curve at 2 K. At higher temperatures there is a more abrupt 
magnetization transition to quasi-saturation by the applica- 
tion of magnetic field of the order of 100 Oe or less. 

small kink develops in the M — H curves but contrary 
to the H || a*,c cases it shifts towards lower values as 
temperature is decreased. On further cooling below Tc 
a narrow, almost rectangular hysteresis opens up around 
zero, indicating a formation of ferromagnetic domains. 
Initial magnetization also shows ferromagnetic charac- 
ter, achieving the saturation value in the virgin curve by 
applying only 100 Oe at 10 K. 

Two characteristic features of a hysteresis loop are re- 
manence and coercivity. Below Tc the remanent moment 
is practically constant with a value /x rem « 0.625/is (per 
one Co ion). On the other hand, the coercive field shows 
a linear dependence on temperature in the log-lin plot, 
as indicated in the inset of Fig. [13] The relation 

H C (T) = H° c ■ e- aT (6) 

has been found to describe well the behavior of many 
nanostructured magnetic systems, like thin magnetic 
filmsi 7 - and amorphous system o 19 ' 20 ! 21 ' 22 . The common 
feature in those systems was the presence of magnetic 
clusters with the well defined anisotropy barrier, where 
jumps of magnetic moments of the clusters over the bar- 
riers are temperature assisted. 

The prominent feature of ferromagnetic order in 
Co7(Te03)4Brg is a very sharp transitions between the 
two saturation states, giving rise to almost rectangular 
hysteresis. Very often, rectangular hysteresis is observed 
in multilayer a 25 ' 26 where magnetic and nonmagnetic lay- 
ers are stacked on top of each other (Co and Pt for exam- 
ple 2 ^). Besides the choice of the constituent materials, 
magnetic properties of multilayers depend strongly on 
thickness of the individual layers, as well as on the growth 
process. A sudden reversal of magnetization occurs when 




T(K) 



FIG. 13: Temperature dependence of the coercive field for 
b-axis field direction. Solid line is guide for the eye only. 



for a critical field one nucleation site for a reversed do- 
main is generated and the avalanche effect is propagated 
throughout the material via strong exchange coupling. 
To the best of our knowledge, Co7(Te03)4Br 6 is the first 
'non-multilayer' material exhibiting the effect of rectan- 
gular hysteresis. In view of its layered structure (stack 
of bc-plane layers), embedding the ferromagnetic compo- 
nent within the layers, the latter interpretation seems at 
least as a consistent possibility. 

Alternatively, the rectangular hysteresis might actu- 
ally be underlined by Stoner-Wohlfarth single-domain 
mode l 23 i 27 . Our single-crystalline Co7(Te03)4Br 6 sam- 
ples are certainly to big to represent a monodomain be- 
low Tc, as one easily verifies from the fact that below Tc 
the net magnetization along 6-axis (as well as along any 
other axis) in Hd c = is M = 0. From the initial (virgin) 
curve (Fig. 1121) one notes however that a quasi-saturation 
is achieved already in a very small applied field render- 
ing the sample practically monodomain in any bigger 
fields. For this reason consideration of M — H hysteresis 
in terms of Stoner-Wohlfarth model makes sense. In the 
latter model the total magnetic energy consists of the two 
terms, the energy of uniform magnetization in the field 
sweeping up and down and the uniaxial anisotropy energy 
of magnetization M as it gets deflected from the pre- 
ferred axis. Total energy minimization results with the 
hysteretic AI—H curves, the shapes of which are precisely 
determined by the direction of applied field with respect 
to magnetization/preferred axis^ 7 .. In the geometry of 
magnetic field aligned with the magnetization, being real- 
ized in this study (Fig. [12]), Stoner-Wohlfarth model gen- 
erates strictly rectangular hysteresis loops characterized 
by coercive field equal to anisotropy field He = 2K/M 
(K is a constant of uniaxial anisotropy). At this field 
value the energy of magnetization in ramping applied 
field just overcomes the anisotropy energy enabling spin 
reversal to take place. 

In the particular case of Coy(Te03)4Br6 it is interest- 
ing to note that the effective preferred axis is oriented 
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approximately along c-axis, not along the magnetiza- 
tion axis (&-axis in our case) as is the case in the origi- 
nal formulation of Stoner-Wohlfarth model. In the case 
of Co7(Te03)4Br6 one has to bear in mind that ferro- 
magnetic moment is just a component of canted three- 
dimensional order thus representing, as pointed out in 
this article's title, a byproduct of global antiferromag- 
netism. It is also interesting to compare values of the 
two critical fields H$f and He at the same tempera- 
tures. H$f is systematically bigger than He for, approx- 
imately, an order of magnitude. This finding is consis- 
tent with elementary understanding of dynamics of mag- 
netically ordered systems: To manipulate ferromagneti- 
cally ordered spins one has to apply field at the order of 
anisotropy field Ha while to do the same with antiferro- 
magnetically ordered spins one has to apply much larger, 
exchange-enhanced field (HaHe) 1 ^ 2 (one bears in mind 
that He ^ Ha)- The latter aspect has recently been 
pointed outi^ in the context of requirements for ultrafast 
spin reorientation technologies. 



IV. DISCUSSION 

In order to discus the results let us first recapitulate 
the main observations presented herewith. 

At Tat, a long-range magnetic ordering of 
Co7(TeC>3)4Br6 takes place, with a temperature 
dependent and incommensurate propagation wave 
vector— kicM- In direct space the latter temperature 
dependence might be associated with global spin re- 
orientations giving rise to ferromagnetic component 
along the &-axis, setting in at Tc- Antiferromagnetic 
response along the two other perpendicular axes (a*- 
and c-axis) is kept unchanged, however. Upon lowering 
the temperature below Tc, the commensurate (CM) 
structure is stabilized and coexists with the ICM 
structure in the short temperature interval. The anti- 
ferromagnetic backbone of the CM structure provides 
a natural explanation for the spin flop (or spin flip) 
transitions and rectangular-shaped hysteresis loops 
observed in sweeping the dc field along a*, c-axis, and b 
axis, respectively. When the first nucleation center for 
the reversed domain is created, the backbone structure 
becomes unstable relative to the 180° rotation of all 
the moments, reversing the total magnetization in an 
extremely narrow field interval (~ 100 Oe at 20 K and 
- 400 Oe at 5 K). 

In an attempt to interpret these observation one has to 
point out that low C2/c symmetry and 4 different (dis- 
torted) environments for magnetic ions make the mod- 
eling for Co7(Te03)4Br6 extremely difficult. Instead of 
attempting the latter here we just elaborate magnetic in- 
teractions and mechanism found responsible for complex 
magnetic ordering of Co7(Te03)4Br 6 . Most obviously, 
there are at least two equally important interactions rul- 
ing magnetism of this compound: exchange interaction 
and single ion anisotropy energy. 



Exchange interactions provide a necessary framework 
for magnetic ordering to set in. As the order established 
at Tjy is incommensurate the exchange interactions Jij 
between magnetic moments of Co 2+ ions necessarily ex- 
tends beyond nearest neighbors, differing in size as well 
in sign. Single ion anisotropy plays pronounced role in 
ruling anisotropic paramagnetic susceptibility. In the or- 
dered phase it provides dominant, or at least sizable, con- 
tribution to macroscopic magnetocrystalline anisotropy, 
K. (Other possible contributions to K rely on magnetic- 
dipolar anisotropy and exchange-interaction anisotropy, 
see, e.g. Ref. [^). In the ordered phase K underlays the 
phenomenon of spin fl(o)ip, introducing the temperature 
through temperature dependence of magnetization^. 

Now, just on ground of specific exchange interaction 
network Jy compatible with the incommensurate order 
and the explicit temperature dependence of magnetocrys- 
talline energy one is able to interpret, at least qualita- 
tively, several important features of Co7(Te03)4Br6 in 
the ordered state, like the temperature dependence of 
kicM and a sudden appearance of ferromagnetism at Tc ■ 

Depending on whether the exchange integrals Jij are 
considered as isotropic or anisotropic there are two pos- 
sible scenarios. In the first, one notes a striking similar- 
ity of magnetic ordering patterns of Co7(Te03)4Br6 and 
rare-earth metals Dysprosium (Dy) and Terbium (Tb). 
The latter elements acquire incommensurate-helical or- 
der at their particular T/v featuring temperature depen- 
dent kicM and ferromagnetic transition, taking place 
a few degrees below Tjy. By keeping aside itinerant- 
magnetism peculiarities, generally important for Dy and 
Tb, the ordering phenomenology itself can convincingly 
be interpreted^ on basis of an interplay between the 
isotropic exchange integrals and temperature dependent 
magnetocrystalline energy. Closer examination of the 
driving force responsible for ferromagnetism of Dy and 
Tb showed however that magnetostriction plays perhaps 
a more direct role than the temperature dependence of 
single ion anisotropy 3 ^. Whether magnetostriction sets 
in at T c in Co7(Te03)4Br6 as well is not fully resolved as 
yet. A high-resolution diffraction study is needed to de- 
tect magnetostriction-related structural changes in vicin- 
ity of T c . 

In the second scenario the exchange interaction is con- 
sidered as being anisotropic, either due to directional 
dependence of exchange integrals or due to antisym- 
metric form of the related spin-spin operator in the in- 
teraction Hamiltonian. In the former case the scalars 
Jij are replaced by the appropriate tensor. Compe- 
tition of anisotropic exchange with uniaxial single ion 
anisotropy as analyzed in mean- field approximation 39 , 
gives rise to different types of possible spin orders min- 
imizing the Gibbs free energy, allowing for spin reorien- 
tations below corresponding temperatures^ -. In the lat- 
ter case anisotropic exchange interaction is manifested 
as an antisymmetric (Dzyaloshinsky-Morya, DM) spin- 
spin interaction. The DM interaction 3 ^, is directly re- 
sponsible for numerous cases of systems revealing weak 
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ferromagnetism in canted antiferromagnet a 1 1 ' 23 . As fer- 
romagnetism in Co7(Tc03)4Br 6 is certainly of canted 
type (Fig. |4|), emerging abruptly from global antiferro- 
magnetism, it is plausible to assume relevance of DM 
interaction, for Coy(Te03)4Br6 as well. With this re- 
spect the mechanism ruling the spin reorientation in 
Co7(TeC>3)4Br6 could be closely related to spontaneous 
spin- flip (Morin transition), taking place deep in anti- 
ferromagnetic phase of hematite, a — Fe.-2.Oz, being ac- 
companied by a DM-based ferromagnetic component. A 
driving force for the Morin transition is identified in com- 
petition of single-ion anisotropy with long-range dipolar 
anisotropy terrr£. 

As pointed out by Becker and coworkers 8 , 
Co7(Te03)4Br6 is composed of chains running along the 
b axis. The chains contain Co(2) and Co(3) ions and 
are interconnected with Co(l) ions to form the be layers. 
The layers are linked through the Co(l) - Co(4) - Co(l) 
connection. As indicated by AOM calculations, Co(l) 
and Co(4) ions are well described in the framework of 
single-ion anisotropy even in the low temperature limit. 
On the other hand, Co(2) and Co(3) ions seem to be 
substantially influenced by the exchange interactions 
through the ligands, indicating a good connection along 
the chain direction and rather weak perpendicular to 
them. The preliminary inelastic neutron scattering 
experiments along the c direction indicate the presence 
of the dispersionless mode around 4 meV, pointing to a 
weak coupling along the c direction, corroborating the 
AOM results. 



V. CONCLUSIONS 

In its ground state Co7(Te03)4Br 6 is a three- 
dimensional canted magnetically ordered system reveal- 
ing antifcrromagnetically compensated sublattices in all, 
but the 6-axis direction. The low temperature mag- 
netic structure is stabilized through temperature depen- 
dent incommensurate wave vector accompanied by an 
abrupt emergence of ferromagnetic component along b- 
axis. Magnetic susceptibility studies demonstrate ex- 
treme anisotropy characterizing Co7(Te03)4Brg. Al- 
though the ferromagnetic component establishes strictly 
along the 6-axis, the effective preferred axis, imposed by 



a sizable single ion anisotropy term of octahedrally co- 
ordinated Co 2+ ions, is directed approximately along c- 
axis. Hidden in a three-dimensionally canted spin ar- 
rangement the latter axis represents, as clearly shown in 
susceptibility and M — H studies, a real 'backbone' of 
antiferromagnetic order. Accordingly, magnetic field (or 
its component) has to be applied along this axis to induce 
entirely reversible spin flop transition. Closer examina- 
tion of M — H scans taken along c- and a*- axis shows 
that the effective preferred axis is actually declined from 
the c-axis for an angle of approximately 15°. From the 
extreme sharpness of the transition it is concluded that 
a phenomenon of spin flip (instead of a spin flop) bet- 
ter describes the observations, classifying Co7(Te03)4Br 6 
into the category of metamagnets. Ferromagnetic re- 
sponse, restricted to the direction of 6-axis only, has been 
related to the phenomenology of multilayers and/or of 
the Stoner-Wohlfarth model due to strikingly rectangular 
hysteresis loops. Co7(Te03)4Br6 obviously represents a 
remarkable magnetic system manifesting competition of 
various magnetic interactions. In the competition there 
is primarily a complex network of exchange interactions 
and a single ion anisotropy energy. Most probably the 
single ion anisotropy represents just one possible compo- 
nent of a more complex magnetocrystalline anisotropy, 
relevant for energy balance in the ordered phase, and in 
the article other alternatives are discussed along the lines 
of knowledge accumulated during the decades. 
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